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Context. We report new constraints on the galaxy luminosity function at z ~9 based on observations carried out with ESO/VLT 
FORS2, HAWK-I and X-Shooter around the lensing cluster A2667, as part of our project aimed at selecting z ~ 7-10 candidates 
i-C ' accessible to spectroscopy. Only one J- dropout source was selected in this field fulfilling the color and magnitude criteria. This 

source was recently confirmed as a mid-z interloper based on X-Shooter spectroscopy. 

Aims. The depth and the area covered by our survey are well suited to set strong constraints on the bright-end of the galaxy luminosity 
function and hence on the star formation history at very high redshift. 

Methods. The non-detection of reliable J- dropout sources over the ~36arcmin 2 field of view towards A2667 was used to carefully 
determine the lens-corrected effective volume and the corresponding upper-limit on the density of sources. 

Results. The strongest limit is obtained for ®(Mi 50 o=-21.4±0.50)<6.70xlO~ 6 Mpc" 3 mag" 1 at z ~9. A maximum-likelihood fit of 
the luminosity function using all available data points including the present new result yields M* >-19.7 with fixed a=-1.74 and 
O* = 1.10xl0~ 3 Mpc~ 3 . The corresponding star formation rate density should be Psfr <5.97xl0~ 3 M o yr 1 Mpc 3 at z ~9. These results 
are in good agreement with the most recent estimates already published in this range of redshift and for this luminosity domain, 
f — ■ Conclusions. This new result confirms the decrease in the density of luminous galaxies at very high-redshift, hence providing strong 

QO ' constraints for the design of future surveys aiming to explore the very high-redshift Universe. 
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, 1 . Introduction observing program combining these two approches aimed at se 

lecting z ~ 7-10 candidates accessible to spectroscopy, namely 
Understanding the formation and evolution of the first galax- the WIRCAM ultra Deep Survey (WUDS) at CFHT on the 
les is one of the major topics of modern astronomy. Since a CFHTLS-D3 Groth Strip field (Pello et al., to be submitted), and 
dozen years, many projects and instruments have been devel- a lensing suryey dLanorteetaLll20lTl hereafter Lll), which is 
opped to push the boundaries of the observed Univer se, such cen t ra i to this letter 
as the Hubble Ultra Deep Field ( Beckwith etal] 120031) for ex- 
ample. The role played by the first luminous objects in the The g alax Y Luminosity Function (LF) is relatively well con- 
Universe during the cosmic reionization is still unclear. One way strained up to z ~6. Beyond this limit however the contamination 
to explore the first epochs is to constrain the number density by low-z interlopers becomes more severe because a majority 
of galaxies along the cosmic time, and thus establish the Star of samples are only supported by photometric considerations, 
Formation Histoiy (SFH) which represents the evolution of the and therefore the results a t these redshifts cou ld be seriously af- 
Star Formation Rate density (SFRd) over the cosmic time. Two fected by contamination dBoone et alj|201 1|). The only way to 
complementary approaches could be adopted to further con- 8 et rid ° f mis problem is to spectroscopically confirm, if not 
strain the luminosity function : large blank fields to select the alL at leas t a substantial fraction of sources in a given sample, 
brightest sources, and lensing clusters used as gravitational tele- a performance beyond the capabilities of present-day facilities, 
scopes to select the faintest ones. Our group has conducted an O nl y the brightest samples, including highly magnified sources, 

could be presently confirmed. Fortunately, the intermediate-to- 

* Based on observations collected at The European Southern bright accessible domain (M uv —20.5 to -22) is also the most 

Observatory, Paranal, Chile, as part of the ESO 082.A-0163 and 087.A- sensitive region of the LF, where most of the evolution is cur- 

0118 rently observed as a function of redshift (e.g. Bouwens et al. 
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2010). However, at z ~7-8 the LF started to be set with recent 
observations from ground based or space observatories (such a s 
lOesch et alj|2012t iMcLure et al.ll201 ll ICastellano et al.l 120 1 Oh . 
Some attempts have been tried up to z ~10 but without spectro- 
scopic confirmation dBouwens et al.ll201 lal) . 

Deep images of the lensing cluster A2667 were obtained 
with ESO/VLT FORS 2 and HAWK-I in /, z, Y, /, H and 
Ks bands. All details regarding these observations and data 
reduction can be found in Lll. IRAC data between 3.6 and 
8yum, and MIPS 24//m were also inclu ded in the analysis 
when avail able (Spitzer Spac e Telescope, IWerner et al.l I2004L 
iFazioetalJ 12004 iRieke et alj|2004l) . Based on these data, 10 
photometric candidates at z ~7-10 were selected based on the 
usual dropout technique. Among them, only one candidate, 
Jl, is consistent with the z~9 /-drop selection criteria, with 
observed rritf=25.21±0.08 and a modest magnification factor 
yU= 1 .3. Despite an excellent fit of the photometric SED at high-z, 
the reliability of this source was already considered as dubious 
by Ll 1 given the detection of a faint compact counterpart on 
the HST F850LP/ACS image, with z 850 =27.39+0.18. Jl was 
recently confirmed as a z=2.082 interloper based on the detec- 
tion of 5 emission lines ( [0111^5007,4959, Ha, Hff and Lyq) 
with ESO/VLT X-Shooter spectroscopy (see lHaves et al.l [2012. 
submitted, for more details). 

In this letter we report on new constraints on the LF at z ~9 
based on the non-detection of realiable /-dropout candidates 
over the field of view towards the lensing cluster A2667. The 
largest magnification regimes cannot be addressed with a single 
cluster because the effective surfaces/volumes explored in this 
domain are too small to retrieve significant results. Therefore, 
this study is limited to the low to moderate magnification regime, 
where this pilot study provides the most interesting results on 
the LF as compared to previous studies, in particular blank-field 
surveys. In Sect. [2] we introduce the method used to constrain 
the bright-end of the LF at z ~9, in particular to determine the 
lens-corrected effective volume and the corresponding upper- 
limit on the density of sources. The constraints derived on the 
SFH at z ~9 are presented in Sect. [3] A brief discussion and 
conclusions are given in Sect. |4] Throughout this paper, a con- 
cordance cosmology is adopted, with £2a = 0.7, Q. m = 0.3 and 
Hq = 7 km s~ l Mpc~ l . A ll magnitudes are given in the AB 
system dOke & Gunnlll983l) . 



2. The luminosity function at z ~9 

The non-detection of realiable /- dropout candidates over the 
field of view toward A2667 is used in this Section to carefully 
determine the lens-corrected effective volume at z ~9 and the 
corresponding upper-limits on the density of sources. In the fol- 
lowing, we only consider the field of view corresponding to the 
surface covered with more than 75% of the total exposure time 
in all filters. A mask was also applied to remove all the noisy re- 
gions within and around galaxies from the subsequent analysis. 
The region masked is only ~ 12% of the total area, but it reaches a 
maximum of ~30% in the cluster core, within the larcmin 2 cen- 
tered on the cD galaxy. The final effective area is ~36arcmin 2 . 

The selection criteria adopted by Ll 1 for 8.5 < z < 9.5 can- 
didates require a detection above 5<x level within a 1.3" diam- 
eter aperture in H, i.e. H < 26.22, and the usual color selec- 
tion as follows: / - H > 0.76, H - Ks < 0.5, and / - H > 
1 .3 x (H - Ks) + 0.76. In practice, the depth in the /-band lim- 
its the reliability of the /-dropout selection to H < 26.0, cor- 
responding to a completeness level higher than 90%. Absolute 



the luminosity function at z ~9 



Table 1. Upper limits for the luminosity function at z =8.5-9.5, 
for different magnifications regimes and UV luminosities. 
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magnitudes in the ~1500A domain (M1500) are derived from 
//-band magnitudes assuming a flat spectrum. 

As compared to blank fields, data obtained in gravitational- 
lensing fields need to be corrected for both magnification and 
dilution effects. A magnification factor fi introduces an enhance- 
ment on the observed luminosity, depending on the redshift of 
the source and the location of the image on the sky, without mod- 
ifying the color-selection window. Dilution reduces the effective 
surface covered by the survey, in such a way that a pixel affected 
by a magnification factor fi has an effective area reduced by the 
same factor on the corresponding source plane. We have used the 
lensin g model for A2667 originally obtained by ICovone et al.l 
(2006) to compute the magnification map at z ~9 using the 
public sof tware Lenst oo^j including the MCMC optimization 
method of lJullo et al.l d2007l) . We explicitly compute the effec- 
tive surface probed by the survey through a pixel-to-pixel inte- 
gration of the magnification map, after masking the pixels lying 
in the object's mask, as explained above. The effective surface 
is 32.55 arcmin 2 , corresponding to a covolume of 57602.5 Mpc 3 
between z ~8.5 and 9.5. 

The effective limiting magnitude for the detection of z ~9 
candidates, however, depends on the magnification factor fi. To 
take this effect into account when deriving upper-limits on the 
density of sources, we have considered different magnification 
regimes, with effective surface/covolume decreasing with in- 
creasing yu. Table Q] summarizes these results and provides the 
non-detection constraints obtained on the LF. Areas, volumes 
and limiting magnitudes quoted in this table are effective lens- 
corrected ones. Given the low to moderate magnification regimes 
explored in this study, uncertainties in the effective surveyed 
volumes due to lensing-modeling uncertainties are considered 
negligible (see also Maizy et al. 2010 for details). The strongest 
limit is obtained for (5(M I 5oo=-21.4±0.50)<6.70xlO- f, ]VIpc- 3 at 
68% confidence level assuming Poissonian statistics. Figure Q] 
displays the independent data points on the LF, together with 
prev ious results from th e literature. We also computed points 
from lZheng et al.l d2012l) . for comparison purpose, using MC ap- 
pro achJakirig^Jnto_jccomit the probability distribution in red- 
shift dBolzonella et al.l l2002) and considering two different val- 
ues for the magnification : fi ~15, as reported by the authors, 
and n ~8, as ob tained from the MACS 1149+22 model by 
ISmithetal.ld2009h . 

We have ado pted the Schechter parametrization of the LF 
(ISchechterlll976h : 

0(M 150 o) = (D*^(lO- a4 < M - M *f +1 exp( - lO" ^"^) 

(1) 

to set constrai nts on M* using a y 2 minimization, assuming a 
fixed a =-1.74 (Bo uwens et alj2008l) . Our two independent lim- 
its in O at M1500 =-21.35 and -20.44 have been combined to the 
available upper limits from the HUDF and GOODS fields at the 

1 http : //www . oamp . fr/cosmology/lenstool 
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Figure 1. Luminosity function at z ~9 showing the upper lim its resulting from th e absence of /-dropout sources in our HAWK-I 
survey (red dots including cosmic variance), upper limit from Bouwens et al. (2008) found in the HUDF and GOODS field of view 
(green dot). We also plot points obtained from the z ~9 candidate of Zhen get alJd2012l) considering a magnification factor of fi ~ 15, 
and the amplification factor /j ~8 obta ined using the model b y lSmith et alj d2009|) for MACS 1 149+22 (error bars take into account 
the cosmic variance computed from lTrenti & Stiavelli] (120081) ^ and points from lLorenzoni et"aTI (1201 ll) at z ~8-9. The last data set 
is only displayed for comparison purposes; it was not used for the computation of the LF at z ~9 given the difference in the redshift 
domain. The dark line displays the Schechter function using M*=-19.7, ar=-1.74 and <t>* = 1.10xl0~ 3 Mpc~ 3 mag~'. The small panel 
displays the likelihood contours for the 1, 2 and 3 cr confidence regions for a fixed a--l.lA. 



same redshift domain than in iBouwens et al.l d2008l) . As shown 
in Fig. [TJ given the degeneracy between $>* and M*, and the 
small number of available data points at z ~9, we can only de- 
rive a limit for <£* when assuming M* or vice-versa. It has been 
shown that <t>* seems not evolve significatively from z ~4 to 9, 
wherea s M* clearly decreases over z ~4-8 (e.g. IBouwens et alj 
(1201 lbh l. Assuming this no-evolution, we have used a fixed 
O* =1.10xl0~ 3 Mpc- 3 and we obtain M* >-19.7. This result is 
also consistent with previous findings (IBouwens et alj d2008l) ). 

3. The star formation history up to z ~9 

The SFH could be explored through the evolution of the SFRd 
over the cosmic time. Using the shape of the LF found in the 
previous section, we are now able to constrain the SFRd at z ~9. 
This density is deduced from the U V density produced by Lyman 
Break galaxies and usually defined by : 

PUV = I il500 ( l ) (il500)^1500 (2) 

Jo.05L* =3 

where L* =3 =1.33xl0 41 erg s _I A" 1 and <t>(Li 500 ) is the UV LF. 
The lower integration bound was chosen in order to facilitate 
the comparison with previous studies (see Sect. |4). T he conver- 
sion b etween UV density and SFRd is done using fhe lKennicuttl 
(1998) relation given by : 

p SF/? (M .yr 1 .Mpc- 3 ) = 1.40 10~ 28 xpt/vKerg.s-'.Mpc^.A- 1 ) 

(3) 



As a result, the SFRd limit at z ~9 is found, using M*=- 
19.7, to be p SFR <5.97xlO~ 3 M yr _I Mpc 3 . Fig.|2]shows the 
evolution of the SFRd versus cosmic time including results al- 
ready published over a large range in redshift. For consistency, 
all points in this figure up to z ~1 have been corrected for 
dust extinction (at z >8 we considered that dust attenuation 
is negligible). A parameterization as a function of redshift has 
be en derived based on t he ex pression for dust extinction given 
by ISchiminovich et all d2005l) as a function of the UV slope, 
the evolution of which as a function of redshift is taken from 
ISchiminovich et all (2005) fo r z <3.5 , and ex tended fromz ~3.5 
up to 7 according to IBouwens et ail d201 lbl) . The expressiorfl 
given by ICole et all d200ll) is used to fit the SFH, in order to 
compare the above value of psfr with the trend observed be- 
tween z ~3 and 8. The best-fit parameters of the Cole function 
are found to be a=0.0, b=0.10, c=2.85 and d=4.65. As seen in 
Fig. 12 our result is compatible with the evolution observed be- 
low z ~8, a nd it is also con s istent w ith the SFRd limit published 
atz -10 bv IBouwens et all d201 lal) . 



4. Discussion and Conclusions 

Over the ~36arcmin 2 field of view of our survey around the lens- 
ing cluster A2667, only one /-dropout candidate was retained 
(Lll). This source, already considered as a possible interloper 



2 P* = irfph M .yr-'.Mpc- 3 
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Figu re 2. Star formation history including t he present result at 
IWvder et all d2005l) iReddv & SteideJd2009h . lvan"der Burg et all 
d2010l) . lOesch et al.ld2012l) ). For consistency, all points in this 
figure have been corrected for dust reddening (see text). The red 
line shows the pa rametrization of the corrected SFH following 
ICole et al.l d200ll) (see text) and the dark line displays the evolu- 
tion of uncorrected SFRd (i.e. without dust extinction correction; 
the corresponding points are not shown to avoid confusion). 

by Lll, has been recently confirmed at z=2.082 based on the 
detection of five emission lines with ESO/VLT X-Shooter spec- 
troscopy (see Hayes et al., submitted, for more details). The non- 
detection of /-dropout candidates in o ur survey is consistent 
with previous non-detections reported bv lBouwens et al.l (T2008) 
in the HDF, HUDF and GOODS fields (4.8, 9.1 and 9.3 arcmin 2 
respectively), with m// (5cr, 0.6arcsec aperture) typically rang- 
ing between ~27.6 to 26.6, and also complementary in terms of 
surface covered and effective depth. 

Assuming that there is no reliable z ~9 candidate in our sur- 
vey up to m#=26.0+2.5xlog(jLf), we take benefit from the depth 
and large area covered to set strong constraints on the bright- 
end of the LF and hence on the SFH at very high redshift by a 
careful determination of the lens-corrected effective volume and 
corresponding upper-limits on the density of sources. We have 
considered different magnification regimes, with effective sur- 
face/covolume decreasing with increasing fi. The strongest limit 
is obtained for 0(M 150 o=-21.4+0.50)<6.70xlO" 6 Mpc- 3 mag- 1 
at 68% confidence level assuming Poissonian statistics. 

We combined our two independent LF limits at M1500 =- 
21.35 and -20.44 with the available upper limits for t he HU DF 
and GOODS fields, extracted from IBouwens et al.1 d2008l) . to 
derive the Schechter parameters for the LF. Assuming a fixed 
a =-1.74 and no-evolution in <D* = 1.10xl0~ 3 Mpc~ 3 , we ob- 
tain a constrain on M * >-19.7. This r esult is consistent with 
previo us findings, e.g. IBouwens et al.l d2008l). ICastellano et al] 
d2010h (M*(z -7)^-20.24) and lOesch et al] d2012l) (M*(z. ~8)^- 
19.80). The corresponding SFRd given the previous LF should 
be p SFR <5.97xlO _3 M yr~' Mpc 3 at z ~9. These results are 
in good agreement with the most recent estimates already pub- 
lished in this range of redshift and for this luminosity domain, 
and confirm the decrease in the density of luminous galaxies 
from z ~6 to 9, hence providing strong constraints for the de- 
sign of future surveys aiming to explore the very high-redshift 
Universe. 

Several projects such as CL ASH dPostman et al .1 l20Llh or 
CANDELS (Groginet aD l201 ll) are currently ongoing to con- 



strain the evolution of galaxies beyond z ~6. However, all sam- 
ples highlighted by these projects need to be confirmed by deep 
spectroscopic observations beyond the limits of current spectro- 
graphs. Near-IR spectrographs are also needed to understand 
the nature of (extreme) mid-z interlopers presently found in 
such deep surveys. T he arrival of multi-object near- IR faci lities 
such as EMI R/GTC dBalcells et al.l 120001 IPello et al.1 1201 lh or 
KMOS/VLT dSharples et al.ll2006llSharples et al.ll2010l) will in- 
troduce a substantial progress in this area. Given the results pre- 
sented in this paper, increasing by a factor of ~10 the sample of 
lensing clusters with deep multi-wavelength photometry over a 
large field of view, similar to our survey in A2667, is expected 
to set crucial constraints on the intermediate-to-high luminos- 
ity domain at z~7-10, and substantially improve the selection of 
spectroscopic targets. Such a cluster sample will provide a ref- 
erence frame allowing to connect the results obtained in blank 
fields with those obtained in the strongest magnification regimes, 
as shown in this letter. 
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